Coherent single atom shuttle between two Bose-Einstein condensates 
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We study an atomic quantum dot representing a single hyperfine "impurity" atom which is co- 
herently coupled to two well-separated Bose-Einstein condensates, in the limit when the coupling 
between the dot and the condensates dominates the inter-condensate tunneling coupling. It is 
demonstrated that the quantum dot by itself can induce large-amplitude Josephson-like oscillations 
of the particle imbalance between the condensates, which display a two-frequency behavior. For 
noninteracting condensates, we provide an approximate solution to the coupled nonlinear equations 
of motion which allows us to obtain these two frequencies analytically. 
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When two phase-coherent quantum systems are 
brought close together, but are still separated by a tun- 
neling barrier, particle current oscillations are induced. 
The exploration of the Josephson effect, predicted in 1962 
[H , was for many years limited to superconducting mate- 
rials, until in 1997 quantum oscillations through an array 
of weak links in superfluid 3 He-B were observed 0|. 

Dilute atomic gases, which upon lowering the temper- 
ature to the sub-/iK-range acquire phase coherence, al- 
low to study macroscopic coherence effects in a highly 
controlled manner, with essentially single-atom accu- 
racy. Calculations based on the two-mode description 
0| of two weakly coupled condensates 0, [H] have ex- 
emplified the rich dynamics of particle imbalance oscil- 
lations inherent to Josephson junctions between BECs. 
In particular, upon increasing the particle interaction, 
two weakly-coupled condensates enter the novel collec- 
tive state of macroscopic self-trapping, which is due to 
the self- interaction of the tunnelingparticles [1, [f|, and 
has been confirmed experimentally 0, [§| . 

In the following, we explore the indirect particle ex- 
change between two condensates, mediated by a single 
"impurity" atom coherently coupled to the two conden- 
sates, which is located in a tight trapping potential at 
the position of the barrier between the condensates. The 
conventional tunneling channel we assume to be strongly 
suppressed by raising the barrier between the two con- 
densates. The dynamics of the type of impurity we con- 
sider - an atomic quantum dot (AQD) - when coupled 
by laser transitions to a single infinite superfluid reser- 
voir, was studied in 0. When the AQD is coupled to two 
well-separated condensates, it was demonstrated in [IB] 
that when the dot-condensate coupling is smaller than 
or comparable to the intercondensate tunneling coupling, 
the effect of AQD on the particle oscillations between the 
condensates is negligibly small. The specific question we 
address in the present work is if the single impurity AQD 
can act as a coherent "shuttle" between the essentially iso- 



lated BEC reservoirs when particle transfer between wells 
is possible only via the coupling of the dot to the con- 
densates: We term this the "strong coupling" limit. We 
find unexpected behaviour in this strong coupling limit; 
namely, apart from expected small and rapid oscillations 
of the particle imbalance due to single particles going to 
and fro between dot and condensates, large-amplitude 
Josephson-like oscillations between the condensates, me- 
diated by the dot, occur at a smaller frequency. 

Our model system consists of two condensates with a 
large (single-particle) potential barrier erected between 
them, and a single impurity atom, coupled to the con- 
densates by a two-photon Raman transition. Assuming 
a symmetric double- well trap, the system is described by 
the following Hamiltonian, 



H = U[\V 1 (t)\ 4 + \V 2 (t)\ 4 ]- K [V*(t)V 2 (t) + h.c. 
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The above Hamiltonian is valid within the two-mode 
approximation [3y for the total condensate wave func- 
tion, &(r,t) = <?i(£)<M T ') + * r 2W02(r), where the single- 
particle wave functions cii^ (r) (normalized to unity) de- 
scribe the particles localized in their respective wells, 
and #1,2 (t) are time-dependent amplitudes represent- 
ing the tunneling process. The pseudospin (equivalent 
to the two-level system represented by the dot) is de- 
fined by the Pauli matrix vector &(t) = (a x ,a y , a z ), and 
<t + = \{a x + i<j y ) = <r_ is a spin ladder operator. The 
"on-site" interaction between the particles is given by 

Ui = gfdr\Mr)\*, « = - I dr [£(Wi(r)V0a(r)) + 
< / , i( r )^Bpc( r ) < ^ ) 2( T *)] denotes the positive tunneling cou- 
pling [4j], and S is the detuning from the two-photon 
Raman transition coupling a single atom into the 
dot. The corresponding coupling parameter (trans- 
fer matrix element) is T — MIr J dr <t>\{r)<f>d{r) = 
HFIr J dr <j)2(r)(pd(r), where the spatial wave function of 
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the dot <frd( r ) is normalized to unity, and flu is the Rabi 
frequency of the two-photon Raman transition. Note 
that while the overlap integrals of dot-condensate and 
condensate-condensate wave functions are both small and 
comparable in order of magnitude, the strong coupling 
limit of T/k — > oo can be achieved by sufficiently in- 
creasing the Rabi frequency flu, i.e., well above the BEC 
single-particle energies in the overlap region, which enter 
the tunneling coupling k. 

From the Hamiltonian (HJ , we derive the coupled equa- 
tions of motion for the condensate variables tfq^ and 
the pseudospin vector s = (\l r d(t)\tr\'l r d(t)) , respectively 
s±,s z , of the AQD; here, is the temporal dot wave 

function. The condensate equations read (h = 1) 



K&2 +T.S_, 



id t V 2 = U\^ 2 \ 2 ^2 - K&l + TS-, 

while the equations of motion for the pseudospin are 



(2) 



id t s 



-5s. 



id t s z = -2T(V{ + &l)s- + 2T{$i 



*2)S+. (3) 



We now scale time with T 1 , and introduce the following 
set of dimensionless control parameters 



UN, 



o 
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q = — r = — 
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where iV = Ni(0) + N 2 (0) = | (0) | 2 + |<P 2 (0)| 2 is 
the sum of the initial number of particles in the left 
and right wells, respectively; we also employ the scaling 
&i — * &i/y/No, Decomposing and s± into their real 
and imaginary parts, we then obtain seven equations for 
the coupled motion of the two condensates and the AQD 
pseudospin, which we solve numerically. 

We first concentrate on the strong coupling case r — > 0. 
Our major result is that large amplitude Josephson-like 
oscillations of the particle imbalance n(t) — (iVi(t) — 
A^(t))/iVo, with amplitude n(0), can be induced by the 
quantum dot, which can coherently transfer atoms one 
by one from left to right and vice versa even when con- 
ventional tunneling is completely switched off (FigQ](a), 
black curve). With increasing a, macroscopic self- 
trapping, defined by an average (n(t)) ^ 0, occurs (Fig. 
[T](b), black curve). Thus, there is a critical value a c , de- 
pending on j3, such that for a > a c particle imbalance 
oscillations are self-trapped and for a < a c , Josephson- 
like oscillations of the particle imbalance with (n(t)) = 
occur. We have studied in detail the dependence of a c on 
(3 and present the results for 0- and ir— junctions in the 
phase diagram Fig. [2 We observe that, for small 0, a c 
increases linearly in (3. For large /?, a c ~ 1//3, implying 
that the critical interaction U c oc T 2 /S in this limit. 

We stress that the self-trapping crossover obtained here 
is very different from the well-known one Q, as it occurs 
also if ordinary tunneling is blocked. In the conventional 



self-trapping scenario, the latter fact would imply (triv- 
ially) that the system is always self-trapped, with zero 
oscillation amplitude. Here, by contrast, the coupling to 
the AQD can induce Josephson-like oscillations for suffi- 
ciently small a, cf. Fig.[TJ To further emphasize the dif- 
ference to the conventional scenario, varying the particle 
number, we find that for large (3 the critical a c becomes 
essentially independent of No, while for small (3, a c de- 
creases approximately linearly in Nq with increasing Ao 
[keeping n(0) fixed], which is opposite to what one would 
expect for macroscopic self-trapping driven by the total 
interaction energy. 
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PIG. 1: (Color online) Large amplitude Josephson-like oscil- 
lations around zero particle imbalance, n = 0, for No = 1000, 
/3 = 10, initial n(0) = 0.6 and initial phase difference <j>(0) = tt 
for a — 0.01 in (a) and a — 0.1 in (b) [where crossover to a 
self-trapped state has occured for T = 0]. Black solid curve: 
r = 0, red (thin grey) curve: T = 0.05. In (a), we display 
in addition the noninteracting case (a = 0) with V = in 
green (thick grey). We assume throughout our calculations 
that there is initially one particle in the dot, s z (0) = 1 and 
that the initial particle imbalance n(0) = 0.6; time is in units 
ofT" 1 . 



FIG. 2: (Color online) The critical value of the scaled interac- 
tion a = a c for self-trapping of the effective Josephson oscilla- 
tions, as a function of the scaled detuning /3, for a <j>(0) = 0- 
junction (black dots) and </f>(0) = 7r-junction (red squares), 
with iV = 1000, r = 0, n(0) = 0.6. Within the shaded areas, 
the system shows Josephson-like oscillations, while inside the 
white area, it is self-trapped. 



We observe that the dimensionless control parameter 
for the occurrence of the Josephson-like oscillations is the 
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FIG. 3: (Color online) Particle number dependence of induced 
oscillations for vanishing conventional tunneling rate, r = 0, 
for iV = 100 (black, I), N = 300 (red, II), and iV = 1000 
(green, III). Parameters are /3 = 10, a = 1, The initial phase 
difference is <j>(0) = in (a) and </>(0) = tt in (b); time is again 
in units of T _1 . 



ratio of the dot energy over the self-interaction energy of 
the condensates: (3/a = 5/UN . Given that the typi- 
cal mean-field energy UNq ~10 nK for particle numbers 
A*o ~ 1000, for (3/a = 10 1 • • • 10 4 , the necessary detuning 
is of order 8 — £JVo [Hz], which is experimentally feasible. 



If we further increase a (decrease /3/a), the ampli- 
tude of the self-trapped oscillations becomes very small; 
Fig. [3] (green curves). However, the magnitude of oscil- 
lations again increases when the number of particles Nq 
is decreased. In this regime of smaller particle numbers, 
one can clearly distinguish a two-frequency behavior of 
the dot-induced oscillations (which is not discernible in 
Fig. [J) . In order to understand the origin of these two fre- 
quencies, we provide below an analytical derivation of the 
oscillations in the limit of noninteracting condensates. 

We have also analyzed the case of a finite, but still 
small r (Fig.[H red curves). The small rapid oscillations 
on top of the envelope of the oscillations rapidly vanish 
if one increases T from zero, and thus occur in the strong 
coupling limit only. Furthermore, we observe from Fig.Q] 
that the change in the effective tunneling oscillation fre- 
quency, increasing the tunneling coupling T, strongly de- 
pends on the value of a. If we are in a state of ordinary 
oscillations around zero population imbalance, Fig.QJa), 
the oscillation frequency changes strongly; conversely, in 
the self-trapped state, we have nearly no change in os- 
cillation frequency. Furthermore, from Fig.QJb), we con- 
clude that by a slight increase of T, we may switch the 
system from a self-trapped state to one with Josephson- 
like oscillations. 



For noninteracting condensates [11], a = U = 0, an 
analytical approximation to the coupled equations of mo- 
tion is possible, from which we are able to reproduce the 
numerically established dot-induced oscillations in that 
limit. Defining the new variables ip — £q + #2 and 
ijj = ^1—^2, we find immediately from Eqs. <[2j) and 
([3]), that ip decouples according to idtip = nip. Hence, 



tp(t) = ip(0)e tKt , and we are left with the equations 

idtip = —Kip + 2Ts_, idtS- = —Ss— — Tips z , 
d t s z = -4TIm[^*s_]. (5) 

These yield directly d t (\ip\ 2 + s z ) = 0, so that 

\iP\ 2 + s z =C = \iP{Q)\ 2 + s z (0) (6) 

is a conserved quantity. Next, we write ip = \ip\e lv and 
S- = |s_|e lA and utilize that the variables |-0|,|s_|,s z 
are directly connected to each other via the conservation 
of both \ip\ 2 + s z and 4|s_| 2 + s 2 z . We get 



dt\il>\ 
d t <p 



-2T|s_|sin(^-A), (7) 
K -2T^cos(^-A), (8) 



Up to now we have treated Eqs. ([2j) and (|3|) without any 
approximations. In the following, we consider the limit 
|s_| /\ip\ <C 1, which is natural given that we assume 
the condensates to be sufficiently large for the Gross- 
Pitaevskii description to apply. Furthermore, we assume 
that s z oscillates around an average value s z , so that the 
average value of |^| 2 is C — s z , accordingly. Similarly, 
we denote the average value of dtp by u). Differentiating 
\ip\ 2 twice we obtain 



8T 



:> " ' +4:T 2 s z + 2(d t p- k)(k-S). (10) 
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Neglecting the first term on the right hand side and re- 
placing 1^1 2 and dt<p by their averaged values, we get 

d 2 ^ 2 ^4T 2 (C-\ijj\ 2 )(C-- Sz )+2(u>- K )( K -SM 2 , (11) 

which results in the analytical solutions 

= C-s z -A cos[2TVC{t-t )l 

s z = s z +A cos[2TVC{t-t )}, (12) 

and the self-consistency condition 

- 2T 2 s z = (tD - k)(k- 5). (13) 

The quantities \ip\ 2 and s z therefore oscillate around their 
average values with the frequency uj\ — 2T\fC. For 
the remaining determination of the constants Ao,to,s z 
and the average frequency Q, we use from now on the 
simplifying assumption of ip(0) > 0, i.e. ip(0) = 0, 
and s_(0) being real. In that case, the amplitude is 
given by Aq = s z (0) — s z , whereas to = 0. More- 
over, a good approximation for ui is then given as the 
average of two extrema of d t (p in Eq. according to 

^ = -^-(0)]T(fe^ + ^i^). After a 
lengthy, but straightforward calculation we finally obtain 

LU = K - 21 — 5 , (14) 
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which also directly determines s z via Eq. 1|T3|) . The final 
result can then be written as 



<p(t) = \fc-s z -A cos(2T 'VCt) e ia \ 

${t) = tP(0)e-' Kt , where Aq = s z (0) - s z . (15) 

Using that Refofr^*] = |^i| 2 — 1^2 1 2 , the particle imbal- 
ance oscillates according to 



n(t) = 



Aq COs(wit) 



Re[e* i2t V*(0)]. (16) 



It executes small oscillations of the amplitude at the fre- 
quency u>i — 2Tv / C ? , cf. Fig. [21 while the major oscilla- 
tions of the envelope are given by 



= 2n-2T 



(k- S)s z {0) + 2cjis_(0) 
(k-J) 2 +w 2 



(17) 



When T — > 0, we have simply conventional Josephson os- 
cillations at the frequency 2k. Conversely, in the strong 
coupling limit, and with one atom initially in the dot, 
the envelope frequency fl depends in a simple way on 
P : T = 7p+Ic ' w ith a maximum at j3 = 2\fC. Com- 



paring the frequency f2 with our numerical results shows 
excellent agreement. For the set of parameters used in 
Fig.Q](a), for example, the analytically obtained periods 
are 284 T" 1 and 52 T^ 1 for T = and T = 0.05, respec- 
tively. The first period corresponds exactly to the green 
curve. Comparing the latter with the red curve shows 
that the agreement remains very good even for small in- 
teractions (finite a), provided T is also finite. 



in which a standing laser beam and crossed conductors 
on a microchip create the required trapping architecture 
3, see Fig.Sl We note that all atoms feel the opti- 
cal potential used in this setup. However, the magnetic 
potential acting on one particular hyperfine state domi- 
nates the optical potential at the AQD location, provided 
that the "contact interaction gap" Udd (the energy bar- 
rier for double occupation of the dot) dominates the other 
relevant energy scales, Udd ^> T, UNo. 

In conclusion, we have shown that an atomic quan- 
tum dot can act as a coherent "shuttle" between two es- 
sentially isolated BECs, transferring atoms from left to 
right and vice versa such that Josephson-like oscillations 
are established. Using an analytical approximation in the 
noninteracting limit, we obtained explicit expressions for 
the two frequencies characterizing the dot-induced oscil- 
lations. We have numerically established a phase dia- 
gram for non-self-trapped and self-trapped phases of the 
system, and analyzed experimental feasibility. An exten- 
sion of the present work is to consider arrays of atomic 
quantum dots in optical lattices, where phenomena like 
self-trapping [14J and the influence of the array on the 
Mott insulating state can be studied. 

We thank N. Schopohl, P. Treutlein, and I. Carusotto 
for helpful discussions. URF and AP acknowledge sup- 
port by the Institut Henri Poincare- Centre Emile Borel. 



Crossed Conductors = Magnetic Potential 




FIG. 4: (Color online) Possible experimental setup for an 
atomic quantum dot coupled to BECs on a microchip (light 
blue). A standing laser wave (green) creates deep potential 
wells perpendicular to the chip surface, equivalent to an opti- 
cal potential for all hyperfine states of the atoms. Two crossed 
conductors on the chip (yellow) create a tight magnetic po- 
tential for one hyperfine species (red) causing a large gap Udd 
for double occupation of the dot. An additional laser (not 
shown) couples the impurity atom (red) and the condensates 
(dark blue), resulting in the transfer coupling T. 

We finally discuss a possible experimental realization 
of our theoretical proposal. For that purpose, we sug- 
gest a combined microtrap-standing optical wave setup, 
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